Absence of extended states in a ladder model of DNA 
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We consider a ladder model of DNA for describing carrier transport in a fully coherent regime 
through finite segments. A single orbital is associated to each base, and both interstrand and 
intrastrand overlaps are considered within the nearest-neighbor approximation. Conduction through 
the sugar-phosphate backbone is neglected. We study analytically and numerically the spatial extend 
of the corresponding states by means of the Landauer and Lyapunov exponents. We conclude that 
intrinsic-DNA correlations, arising from the natural base pairing, does not suffice to observe extended 
states, in contrast to previous claims. 
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I. INTRODUCTION 

According to standard theories of disordered systemsji 
all states in low-dimensional systems with uncorrelated 
disorder are spatially localized. Therefore, in a pure 
quantum-mechanical regime, DNA might be insulator 
unless the localization length reaches anomalously large 
values. To explain long range electronic transport found 
experimentally, 2 Caetano and Schulz claimed that intrin- 
sic DNA-correlations, due to the base pairing (A-T and 
C-G), lead to electron delocalization. 3 Furthermore, they 
pointed out that there is a localization-delocalization 
transition (LDT) for certain parameters range. If these 
results were correct, then transverse correlations arising 
intrinsically in DNA could explain long range electronic 
transport. However, we have claimed that this is not 
the case and all states remain localized, thus excluding 
a LDT<4 Therefore, the electrical conduction of DNA at 
low temperature is still an open question. 

In this paper we provide further analytical and nu- 
merical support to our above mentioned claim, aiming 
to understand the role of intrinsic DNA-correlations in 
electronic transport. To this end, we address signatures 
of the spatial extend of the electronic states by means of 
the analysis of the Landauer and Lyapunov coefficients, 
to be defined below. The outline of the paper is as fol- 
lows. In the next section, we introduce the ladder model 
of DNA4 and diagnostic tools we use to elucidate the 
spatial extend of electronic states in the static lattice. 
In Sec. IIIII we discuss the analytical calculation of the 
Landauer exponent and show that this exponent never 
vanishes in the thermodynamics limit for any value of 
the system parameters. From this result we conclude 
that extended states never arise in the model. We then 
proceed to Sec. IIVI in which we numericaly calculate the 
Lyapunov exponent for finite samples. We discuss in de- 
tail its dependence on the model parameters, especially 
inter- and intrastrand hoppings. We provide evidences 
that the localization length is only of the order of very 
few turns of the double helix for realistic values of the 
model parameters. Therefore, this shows that intrinsic 
DNA-correlations alone cannot explain long range elec- 



tronic transport, found in long DNA molecules^ Finally, 
Sec. IVl concludes the paper. 



II. MODEL AND DIAGNOSTIC TOOLS 

Our analysis proceeds as follows. We consider a ladder 
model of DNA in a fully coherent regime and assign a sin- 
gle orbital to each base. Conduction through the sugar- 
phosphate backbone is neglected hereafter. Both inter- 
strand and intrastrand overlaps are considered within the 
nearest-neighbor approximation. We assume that the 
hopping does not depend on the base and therefore only 
two values are considered, namely interstrand (tj_) and 
intrastrand (ty) hoppings. Figure ^ shows a schematic 
view of a fragment of this ladder model. 
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FIG. 1: Schematic view of a fragment of the ladder model for 
DNA molecules, excluding the sugar-phosphate backbone. A 
single orbital is associated to each base, with a corresponding 
energy e„, n being A, T, C or G. The sequence of the basis 
of one of the strands is random, while the sequence of the 
other strand results from the base pairing A-T and C-G. Only 
two different hoppings are considered, namely interstrand (t±) 
and intrastrand (tii) hoppings. 

Four different values of the energy sites (ea, £t> £c, 
and eq) are randomly assigned in one of the strands, 
with the same probability, while the sites of the second 
strand are set to follow the DNA pairing (A-T and C-G). 
Hereafter we will restrict ourselves to the following val- 



2 



ues of the site energies, taken from Ref. 0, sa = 8.24 eV, 
e T = 9.14 eV, e c = 8.87 eV, e G = 7.75 eV. The same site 
energy values were taken in Ref. 0- As a consequence, 
only three parameters remain in the model, namely tj_, 
in and N, the number of base pairs. We will show below 
that the spatial extend of the states strongly depends on 
ij_ and iii, but never goes to infinity in the thermody- 
namics limit (N — > oo ). 

Once the model has been established, we can write 
down the equation for the amplitudes at different bases. 

Let us denote these amplitudes as ip^ , where a = 1 , 2 
runs over the two strands and n = 1,2, ... ,N denotes 
the position of the bases at each strand. According to the 
model introduced above, the equations for the amplitudes 
are readily found to be 

E*2> = e^g* + t„ (V«r + &) + , (la) 

= 4 2 Vi 2) + *|| + ^i) + t ± ^ . (lb) 

Here takes one of the four values of the site energies, 
according to the constraints presented above. 

The equation for the amplitudes can be cast in a com- 
pact form by using 4x4 transfer matrices. To this end, 
let us introduce the 4-vector 



(tfKtfW&^i)* (2a) 



where the superscript t indicates the transpose. Defining 
the following 2x2 matrix 



M. 



(n) _ 



(2b) 



we arrive at the transfer-matrix equation <f>„+i = T n #„, 
with 



" 1 /, o 2 



(2c) 



where 2 2 and 2 are the unity and null 2x2 matrices, re- 
spectively. One can easily find out, that the 4x4 transfer 
matrix T n satisfies the condition T^JT n — J with 



J = 



/O -i 

-i 

1 
\0 i 



(3) 



the localization length of the states, assuming exponen- 
tial localization^ Here i,j G {A, T, G,C} and length is 
measured in units of the base separation along a sin- 
gle strand (3.1 nm). Due to the self-averaging property 
they can be calculated by taking the product of ran- 
dom transfer matrices over a long system. Similarly, in 
the Landauer exponent jLa(N) = M ||IIn=JV T n,ij\\ ) 1/N 
(hereafter (...) denotes ensemble averages) is twice the 
largest Lyapunov exponent near the critical region in one- 
dimensional system*^ and can be calculated a nalytic ally 
following the technique developed in Refs. l9l llOlllL In 
quasi-one dimensional chains, as in the ladder system un- 
der consideration, both exponents again exhibit the same 
critical behavior (the critical indices are the same), but 
their ratio can be different from 2. 



III. LANDAUER EXPONENT 

In Refs. l9l llCilll it was shown that the Landauer resis- 
tance and the corresponding exponent can be calculated 
exactly. To this end, the direct product M. N ® M^ N of 
the fundamental representations of transfer matrices of 
the SU(1, 1) group is reduced to the adjoint one. We 
apply this technique here for the group SU "(2,2). In or- 
der to calculate this direct product exactly, we use the 
known representation of the permutation operator via 
generators r p (jjl — 1, . . . , 15) of the sl(4) algebra as 
P = (1/4) (I ® I + ® t^). Thus the matrix elements 
satisfy 



Sja 2 X@2 



(4) 



where we assume summation in the repeated indices jU. 
Among of generators r M there is one, which coincides 
with the metric J defined in J3J. We denote the corre- 
sponding index [i as J, namely r J = J. 

Multiplying J3J by Tj and TJ from the left and right 
hand sides respectively, one can express the direct prod- 
uct of Tj and Tj via their adjoint representation 



= l{J)%J) P al + j(r"J)^Af (Jr")1 . (5) 
Here the adjoint representation A„ of T n is defined by 
Ar = -Tr(T„T^T+0 , (6) 
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which means that T n belong to the SU(2, 2) group. It is 
to be noticed that only four transfer matrices appear in 
DNA due to the intrinsic pairing (see Fig. [TJ, and they 
will be denoted as T n , AT , T n ^ TA , T n ^ C G, T n ,GC for the 
sake of clarity. 

The electronic properties can be described by the full 
transfer matrix Mn = Yl n =N^ 1 n,ij and the Lyapunov (Mn-M^) 
exponents, 7L y , the eigenvalues of the limiting matrix 
limjv^oo ln(A / lAr.Mj v ) 1 / 2Ar , provide information about 



being an 15 x 15 matrix that depends on the parameters 
of the model at site n of both chains. 

To calculate analytically the Landauer exponent, we 
apply this decomposition to the products of fundamental 
representations, T„'s in A4jv, and after averaging obtain 



\j®j+\{t»j)®{jt») |n<A^ 



(7) 
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It is then straightforward to get the average over four 
equivalent substitutions of the base pairs in the random 
chain 

(AH = \ (Kt + Ka + «cg + Kc) ■ (8) 

The Landauer exponents are the nonnegative eigenvalues 
of i log( A) . The condition of the existence of an extended 
state is equivalent to det |(A) — I\ = 0, and it is a matter 
of simple algebra to prove that this condition is never 
met. 

Therefore, we come to the conclusion that the system 
studied by Caetano and Schulz- cannot support truly ex- 
tended states. Consequently, a LDT is not to be observed 
since all states are spatially localized. 
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IV. LYAPUNOV EXPONENT 

It can be argued that, although the localization length 
is always finite, as we have demonstrated above, it could 
be larger than typical sizes used in transport experi- 
ments, as those carried out by Porath et alrk To quan- 
titatively determine the spatial extend of the electronic 
states, we have numerically calculated the Lyapunov ex- 
ponents for different values of the hoppings t± and tu. 

Figure |21 shows the Landauer and Lyapunov exponents 
for N = 4000, when t x = 0.5 eV and t\\ = 1.0 eV. These 
values of the interstrand and intrastrand hoppings are 
larger than those usually considered in the literaturejiSii^ 
but they were used by Caetano and Schula^ to provide 
support to their claim about the extended nature of the 
states. From Fig. it becomes clear that neither the 
largest Lyapunov exponents nor the Landauer one van- 
ish over the whole energy spectrum. Most important, its 
minimum value is size independent within the numeri- 
cal accuracy, suggesting the occurrence of truly localized 
states. Notice that the minimum value of these expo- 
nents is always much larger than the inverse of the num- 
ber of base pairs (1/N = 0, 00025), indicating that DNA- 
pairing can hardly explain long range charge transport at 
low temperature. From the inverse of the minimum value 
of the second Lyapunov coefficient we can estimate that 
the localization length is of the order of 80 base pairs (i.e. 
roughly eight turns of the double helix), therefore being 
much smaller than typical sizes used in experiments^ 

To elucidate the effects of the base pairing on the lo- 
calization length, we have also considered the artificial 
case of ladder models without pairing. In that case, 
both strands are completely random, allowing for a larger 
number of possible pairs (e.g. AC or AA). Therefore, the 
system becomes much more disordered and one could 
naively expect a dramatic decrease of the localization 
length, as compared to the system with base pairing. 
Figure [3] indicates that this is not the case. The in- 
verse of the second Lyapunov exponent remains almost 
unchanged over a large region of the energy spectrum, 
except close to the two resonances at about 6.4 eV and 



FIG. 2: 2/3 of Landauer exponent 7L a (solid line) and the 
largest Lyapunov exponent 7L y (white circles) , as a function of 
energy, for t± = 0.5 eV and iy = 1.0 eV. The second, smaller 
Lyapunov exponent (grey circles) is also shown. 



10.6 eV. At resonances the localization length is reduced 
by a factor 2.5 at most when the pairing constraint is 
relaxed. In any event, resonances still appear so they 
cannot be associated to base pairing. 
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FIG. 3: Inverse of the second Lyapunov coefficient for N = 
4000 t± — 0.5 eV and tu = 1.0 eV, when the base pairing is 
present (solid line) and absent (dashed line). 

We claimed that the spatial extend of states strongly 
depends on the hopping parameters^ and those consid- 
ered in Ref. seems to be larger as compared to those 
values widely admitted in the literature^Siii Higher hop- 
pings lead to a less- effective disorder and higher localiza- 
tion lengths are to be expected. We have calculated the 
inverse of the second Lyapunov exponent for more realis- 
tic values of the hopping parameter and checked that this 
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claim is indeed correct. For instance, for t± = 0.05 eV 
and t|| = 0.5 eV (see Ref. El) the localization length at 
the center of the band is reduced by a factor 5 as com- 
pared to the case shown in Fig. ^ while an even larger 
decrease is noticed at resonances. Therefore, we come 
to the conclusion that hopping is more important than 
base-pairing in this ladder model. 



regime. In this model, a single orbital is assigned to each 
base, and the sugar-phosphate channels are excluded. 
Therefore, it is assumed that electronic conduction takes 
place through orbital overlap at the bases. The sequence 
of bases of one of the strands is assumed to be totally 
random, while the sequence of the other strand results 
from the base pairing A-T and C-G. 
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FIG. 4: Inverse of the second Lyapunov coefficient for N = 
4000 and two sets of hopping parameters, indicated on the 
legend box. Notice the scaling factor indicated on the lower 
curve. 



We demonstrated analytically and numerically that 
due to the randomness of the sequence, the states are 
always localized and a LDT cannot take place, contrary 
to what claimed in Ref. H In particular, we observed 
that base pairing has negligible effects on the localiza- 
tion length except close to two resonant energies, lo- 
cated at about 6.4 eV and 10.6 eV for t± = 0.5 eV and 
tii = 1.0 eV. At these particular energies the localization 
length is smaller when the pairing constraint is relaxed. 
Most important, even in the case of base pairing, the 
larger localization length is much smaller than typical 
sizes of samples used in various experiments. Therefore, 
we come to the conclusion that base pairing alone is un- 
able to explain electronic transport at low temperature 
in DNA. 
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